Spin response of a normal Fermi liquid with noncentral interactions 
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We consider the spin response of a normal Fermi liquid with noncentral interactions under con- 
ditions intermediate between the collisionless and hydrodynamic regimes. This problem is of im- 
portance for calculations of neutrino properties in dense matter. By expressing the deviation of 
the quasiparticle distribution function from equilibrium in terms of eigenfunctions of the transport 
equation under the combined influence of collisions and an external field, we derive a closed ex- 
pression for the spin-density-spin-density response function and compare its predictions with that 
of a relaxation time approximation. Our results indicate that the relaxation time approximation is 
reliable for neutrino properties under astrophysically relevant conditions. 

PACS numbers: 21.65.-f. 26.50,+x, 33.25.+k, 67.30.em 



I. INTRODUCTION 

Neutrino scattering, production, and annihilation rates 
are key physical quantities in understanding stellar col- 
lapse, neutron star formation, and supernovae, and they 
are determined by the spin, isospin and density responses 
of nuclear matter [l[ . The challenge is to include the ef- 
fects of collisions between nucleons and mean-field effects 
in a consistent way. This problem has been approached 
in the framework of Landau's theory of normal Fermi 
liquids @ but, in the detailed calculations performed to 
date, the collision integral has been replaced by a simple 
expression with a single relaxation time [3| . 

In this paper we consider the spin response at long 
wavelengths. This problem is well understood in two 
limiting cases. The first is the collisionless regime, where 
the collision rate is small compared with the frequency ui 
of the field, and collisional effects may be included per- 
turbatively fs}. When collisions cannot be taken into 
account perturbatively, they give rise to the Landau- 
Pomeranchuk-Migdal effect, which in the language of 
perturbation theory is due to energy denominators ac- 
quiring an imaginary part due to collisional damping of 
excitations [J|. 

The second regime that is well understood is the hy- 
drodynamic one, where the collision rate is large com- 
pared with the frequency. In this case, the transport 
equation for quasiparticles may be solved exactly using 
the methods developed by Jensen, Smith and Wilkins [5[ 
and by Brooker and Sykes @. In astrophysical applica- 
tions, conditions intermediate between these two limiting 
cases arise. Related problems for the spin-independent 
part of the response, corresponding to distortions of 
the quasiparticle distribution in momentum space with 
an angular dependence proportional to a spherical har- 
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monic, have previously been attacked using variational 
methods. These were developed in the context of the 
frequency-dependent electrical conductivity and magne- 
toresistance of a normal metal Q, and provide upper 
and lower bounds for the imaginary part of the response. 
These methods were subsequently exploited to obtain up- 
per and lower bounds on the response of the quadrupolar 
distortion of the quasiparticle distribution, which is a key 
feature in calculations of the attenuation of sound in the 
regime intermediate between the first and zero sound [8[ . 

Exact expressions for relaxational modes of a Fermi 
liquid have previously been derived in the regime inter- 
mediate between the collisionless and hydrodynamic ones 
in Ref. @. In that paper, eigenfunctions of the Fourier 
transform (with respect to time) of the Boltzmann trans- 
port equation for quasiparticles were found. 

In the present paper we extend the work on exact so- 
lutions to include an external field, in contrast to most 
earlier work with exact solutions which considered trans- 
port coefficients in the hydrodynamic limit, where the 
effects of the external field drop out [1, @], or involved 
the spectra of relaxational modes and collective modes 
in the absence of an external field ||. We use the eigen- 
functions obtained in Ref. [9( to derive a closed expression 
for the spin-density-spin-density response valid generally. 
We compare our results with those of a single relaxation 
time approximation, with a relaxation time chosen to re- 
produce the correct behavior in the collisionless regime, 
and find that for astrophysically relevant conditions, the 
response function given by the relaxation time approxi- 
mation differs from the exact result by less than 10%. 



II. BASIC FORMALISM 

We calculate the response of a normal Fermi liquid to 
a magnetic field in the z-direction that is spatially ho- 
mogeneous and varies in time as e~ 1QJt . We assume that 
the quasiparticle density matrix and the quasiparticle en- 
ergy may be taken to be diagonal in a z . In doing so, 
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we neglect the fact that noncentral components of the 
interaction, which are very important in systems of nu- 
cleons, can create nonzero x- and y-components of the 
spin-dependent part of the quasiparticle density matrix. 
In addition, we neglect tensor components of the mag- 
netic moment. In the presence of a uniform magnetic 
field, these would give rise to a quadrupolar distortion of 
the quasiparticle distribution as a function of direction in 
momentum space [Toj . Therefore we write the change in 
the energy Se pa of a quasiparticle with momentum p and 
z-component of the spin (o = ±1) due to the application 
of the external field U, as 



from interactions with other quasiparticles, 



-pa 



-pa 



Ss 



pa i 



(4) 



where e pcr denotes the quasiparticle energy in equilibrium 
in the absence of the magnetic field and Se P a includes 
the interaction with the magnetic field, Eq. |T]), and the 
contribution due to quasiparticle interactions fi 2 , 



8ei = at U z 



2 



/12 Sn 2 ■ 



(5) 



S£pa = (?U Z 



(1) 



and denote the quasiparticle distribution function as n pa . 

We start from the quasiparticle kinetic equation, which 
for spatially uniform conditions is given by 



dn pa 
dt 



= I\n T 



(2) 



where /[np/ 5 i] is the collision integral. In normal Fermi 
liquids at low temperatures, the dominant collisions are 
those between pairs of quasiparticles and, in a compact 
notation in which p^ is denoted simply by i, Eq. @ 
takes the form (see Ref. [ll[ for details) 

^ = -m 2nJ2l\ <12|-4|34) \\ 2 (1 - n 3 ) (1 - n 4 ) S e S p 

0t 234 1 

+ (1 - u x )27r^ i|(34|^l|12)| 2 (l - n 2 )n 3 n 4 S £ S p , 

(3) 



234 



where (12|„4|34) is the quasiparticle scattering amplitude, 
5 £ 5 P is shorthand for 5(£i+£2— £3— £ 4 )<5(pi+P2— P3— p 4 ) 
and — ~52a- I dpi / (2tt) 3 . The two quasiparticles in 
the scattering process are of the same species and the 
factor 1/2 avoids double counting of final states when 
momenta are freely summed over. In addition, we as- 
sume that lj <C T, so we do not include u in the energy- 
conserving delta function 

The quasiparticle energies include contributions due to 
the external field and due to "molecular fields" resulting 



Tensor components of the Landau quasiparticle interac- 
tion are generally small in nuclear matter and we there- 
fore neglect them. In this case the response of the quasi- 
particle distribution to a spatially uniform magnetic field 
is isotropic and the interaction energy reduces to gaS, 
where S — J2i CT i <^ n i is the total spin and go (or /g in 
the quantum liquids' literature [H|) is the isotropic com- 
ponent of the spin-dependent interaction between quasi- 
particles. The change in the quasiparticle energy is thus 



fci = ai(U z + g S) . 



(6) 



We now linearize Eq. ([3]) about the equilibrium distri- 
bution function 



1 



3 (e?-/*)/T + X ' 



(7) 



in the standard way [Tl|, Sect. 1.2.4] taking into account 
the variation of the quasiparticle energies in the delta 
function. Here 11 is the chemical potential in equilibrium. 
It is convenient to introduce the quantity 

dn° 

Srii = rii - n°(ei) = m - rq - -Q^ Se i > ( 8 ) 



which represents the difference between the quasiparticle 
distribution function and the Fermi function evaluated 
for the actual quasiparticle energy, not the one in equi- 
librium. Linearizing Eq. Q leads to 



dSm , dn\d{U z +g Q S) 



dt 



del 



dt 



|E 5|<12M34>| a M.(n- 



T ^— ' 2 
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1 2 cosh(a;i/2) t 

x [2cosh(a;i/2)Si + 2cosh(a;2/2)S2-2cosh(a; 3 /2)S3-2cosh(a;4/2)S4] , (9) 



where xi — (e? — A*)/T, so that one has (dn® / 'de^)^ 1 = 
—T [2 cosh(xi/2)] 2 , and S x = S(xi + x 2 — x 3 — X4). In ad- 
dition, we have introduced Hj = 2 cosh(xi/2) Srii, which 
simplifies the subsequent calculations. In writing the 



right-hand side of Eq. ^ in this form, we have used 
the property that for interactions invariant under time 
reversal (pi<7i, P2C2|^|P3C r 3 , p 4 (7 4 } = «-pa - 0-3, -p 4 - 
<J4\A\ — pi 01, — p 2 — ct 2 ))*, an d therefore, since the dis- 
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tribution function is isotropic in momentum space, one transforming with respect to time, Eq. ([9]) may be recast 
may combine the two terms in Eq. ([3]). After Fourier as an inhomogeneous integral equation for S, 



cn(Uz +9qS) 
2Tcosh(xi/2) 



= S 1 + i^^i|<12|^|34)| 
u> 1 * — ' I 1 

234 



n 



i 



2cosh(a; i /2) 



x [2cosh(xi/2) S x + 2cosh(x 2 /2) S 2 - 2cosh(a;3/2) 5 3 - 2 cosh(x 4 /2) S 4 ] 



(10) 



For T <C £f, scattering processes are strongest for quasi- 
particle states in the vicinity of the Fermi surface. As 
a consequence, integrals over the magnitudes of the mo- 
menta may be decoupled from integrals over quasiparticle 
energies, and the scattering amplitudes involved can be 
replaced by their values for all quasiparticle momenta on 
the Fermi surface. In the literature on quantum liquids 
it is customary to write [HI, [l4| 



when ^3 is held fixed, or 



d cos 8-> 



<9£4 



Thus we may write 



Pkv 



for fixed p% . 



(15) 



dp2 = p\ dp2 d cos 8 d(f>2 , 

where 8 is the angle between pi and p 2 and 
azimuthal angle, and 



(11) 



dp 2 dp 3 = 



m* i kF 
P 



de% de 3 de^ d cos 8 d<p e?0 2 , (16) 



dp3 = p 2 3 dp 3 d cos 8 3 



(12) 



where 83 is the angle between p 3 and the total momen- 
tum P = pi + P2 of the incoming quasiparticles and 
<t> is the angle between the plane containing pi and p 2 
and that containing P3 and p 4 . The integral over 83 
can be converted into one over £4 by using the fact that 
P4 = P p3 , and therefore 



9 
Pi 



p2 



2 

■Pa 



2Pp 3 COS 6*3 



from which is follows that 



p 4 dp4 = — Pp 3 d cos 8 3 



(13) 



(14) 



where P = 2k F cos(6>/2). 
Next we express S as 



(U z + g S) 
= o-i £[Xi) 



(17) 



substitute this expression into Eq. (|10j) . multiply the re- 
sult by o 1 , average over 01, and perform the integrals. 
The integral over </> 2 gives 2ir and the term involving Si 
is proportional to the integral K 3 , Eq. (|A5|) . One thus 
obtains the linearized quasiparticle transport equation, 



2cosh(iri/2) 2cjtq 



Here we have introduced the quantities 



C 2 = i-i- 



2uiTn 



1 7T T 2 



tq 16 kpvp 



TrJ(12|.4|34)| 



A. = 



Tr ff |(12|^|34)ha 1 (-a 2 + a3 + a 4 ) 



Tr CT |(12|^|34)| 
where the trace Tr CT denotes a sum over all spins (Ji and 



(19) 
(20) 

(21) 



iA CT f°° dx2dxzdx^5 x 
WT o J-00 4cosh(x 2 /2) cosh(cc4/2) 

I 



the average is defined by 

1 r sin 8d8 



(18) 



(...) 



8?r J cos(6»/2) J 



2tt 



(22) 



and vp = k^/nn* is the Fermi velocity. The factor 
of 1/ (8_7r) ensures that the average is normalized to 
unity [151 ]. The first term on the right-hand side of 
Eq. (TTSJ) includes the contribution from the integral over 
X2, x 3 and X4 given by the function K 3 , Eq. (|A5[) . in Ap- 
pendix El The integrals over x 3 and X4 in Eq. (j 18(1 yield 
the standard form for the kernel, given by the function 
K2, Eq. (| A4|) . However, for the Fourier transformation 
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in the next section, it is simpler if the integral is left in 
the original form. 

Physically, the quantity 2to/tt 2 is the relaxation time 
of a quasiparticle in a particular momentum state at the 
Fermi surface, if deviations of the quasiparticle distri- 
bution from equilibrium in all other states vanish. For 
an interaction with only central terms, the expression, 
Eq. (I20p , reduces to that derived in the context of liquid 
3 He [III. The quantity 1— A CT is a measure of the effective- 
ness of collisions in relaxing the spin of the quasiparticle: 
if spin is conserved in collisions, o\ + 02 = 03 + (74 and 
therefore X a — 1 , while if A CT = 0, the integral term in 
Eq. (fT8|) vanishes and the equation becomes an algebraic 
one, whose physical content is equivalent to a relaxation 
time approximation with a relaxation time that depends 
on the quasiparticle energy, as we will see in Sect. IIVI 

In applications it is sometimes convenient to express 
the scattering amplitudes in terms of the momentum 
transfers k = pi P3 and k' = pi P4. For quasi- 
particles at the Fermi surface these satisfy the condition 



P 2 + k 2 



(23) 



The magnitudes of the momentum transfers are given by 
k = 2fc F sin(6»/2)sin(0/2) and k' = 2k F sin(0/2) cos(0/2) 
and the Jacobian for the transformation is 



9(cos 9, 4>) 



d(k,k') 



kp 



(24) 



and therefore [|[ 



(•••) = 



2i:i 



dk 



2kF dk' fcp 6 (4fcf -k 2 -k' 2 ) 
k F V 4fc F -k 2 -k' 2 



(25) 

When effects of the medium on scattering amplitudes are 
taken into account, the scattering amplitude will also de- 
pend on the total momentum of the incoming quasipar- 
ticles, and the natural variables for scattering of quasi- 
particles at the Fermi surface are P and <j>. In terms of 
these the average, Eq. (|22p. takes on the simple form 



(...) 



2 At 



dP 

2~kZ 



2vr 



(26) 



In the next section, we solve Eq. (fT5|) for £, from which 
one can calculate the total spin S, 



9n? 

(U z +g S)). (27) 



1 

—1 



S = <7i Sni = (7i [ Sni + (Ji -jgy (U z + g S) ) , 



^°' 1 l,2msh(.r l /2) 



0"! 



Since S is proportional to U z + QoS, this equation shows 
that S can be written in the form 



S 



(28) 



where the spin-density-spin-density response function Xa 
is given by 



X<7 



1 + g x a 



(29) 



and X a = —S/(U Z + goS) denotes the response func- 
tion in the absence of mean- field effects (go = 0). Using 
Eq. (pH), it follows from Eq. ((2Tj) that 



X rT = N(0)(l- J 



dx 



_ x 2cosh(x/2) 



(30) 



with N(0) — ~J2id n i/d £ i being the density of states 
at the Fermi surface. 



III. SOLUTION OF TRANSPORT EQUATION 



To solve the linearized transport equation, Eq. (fl8j) . we 
follow the methods developed by Brooker and Sykes 
and express the solution £ in the form 



(31) 



r=0 



where $ r (x) are eigenfunctions of the homogeneous equa- 
tion, 



(xl + 7r 2 C 2 )^ r { Xl ) 

A r f 00 dx2dx 3 dx4 8 X 



2 J _ oa cosh(x2/2) cosh(a;4/2) 



*r(*3). (32) 



The eigenfunctions may be determined by Fourier trans- 
forming. We write 



dx & r (x) e' 



-ikx 



(33) 



By using the representation, Eq. (|A3[) . for the delta func- 
tion, one finds that the Fourier transform of the in- 
tegral term reduces to the product of $ r (fc) and the 
square of the Fourier transform of 1/ cosh(x/2), given 
by 27r/ cosh(7rfc). Thus the eigenfunctions satisfy 



d 2 $ r (k) 
dk 2 



27T 2 A r 

cosh 2 (7rfc) 



$ r (fc) = 0, (34) 



which has the same form as the Schrodinger equation 
for a particle moving in a one-dimensional potential ~ 
sech 2 (7rfc) with a (generally complex) energy ~ ^ 2 
This is analogous to the Sturmian method in quantum- 
mechanical problems, in which one uses a basis of states 
with the same energy, but with different strengths of the 
potential p^ |. 

In this paper we are concerned with disturbances that 
are even in x, and therefore the relevant eigenvalues are 



A r = -(2r + C)(2r + l + C) ! 



(35) 
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where r is positive integer or zero, and the corresponding 
solutions [9|, p. 85] are given by 

* r (fc) = 2 C T(l + C) P a "£ c (tanh(Trfc)) , (36) 

- r^cVl] [ COsh ^)]~ C Ci +1/2 (tanh (7 r fc )), 

(37) 

= 2 F X (-2r, 2r + 1 + 2C, 1 + C, [1 - tanh(7rfc)]/2) 
x [cosh(7r/c)]- c , (38) 

where C™ are Gegenbauer polynomials and %F\ is the 
hypergeometric function. The corresponding eigenfunc- 
tions as a function of energy, $ r (x), may be determined 
by Fourier transformation, and some of their properties 
are described in Appendix [XJ 

On substituting Eq. (|3Tj) into Eq. (|18p . multiplying by 
Q r (x\) and integrating over x\ one finds 



-\2lotc) 



U r (A r - A(j ) 



(39) 



where, following the notation of Ref. @, we have defined 

dk * r (fc) _ r(c + 1) r(i + r + |C) r(J + r) 



^ cosh(Trfc) Tr T(l + r + |C) r(l + r + C) 



and the normalizing integral is given by 



(40) 



U r = 27T 



dk [$ r (k)Y 



2 2 < +1 (2r)![r(C + l)] s 



cosh 2 (7rfc) (2r + C+ ^)T(2r + 2C+ 1) 

(41) 

The integral in Eq. (f30|) therefore leads to 



2cosh(a;/2) 
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-VjOTq 



E 



i 



(42) 



and the response function X a in the absence of mean-field 
effects is given by 



N(0) 



= l + i^ V 

7T 2 ^ 



(2r + C + |) T(2r + 2C + 1) / r(± + r + k) r(± + r) 



2 2C+i( 2r )!(A r -A (T ) 



r(i 



|C)r(i 



•0 



l + i 



2^r CT (l - A(j) ^ (2r + C + |) T(2r + 2( + 1) { r(± + r + ±£) r(± + r) 



E 

r=0 



2 2 C+i(2r)!(A r -A ff ) 
I 



r(i 



k)r(i 







(43) 



The sum in Eq. (|43[) may be expressed in terms of the gen- 
eralized hypergeometric function 5F4 but in applications 
it offers no advantages compared with a direct evaluation 
of the sum. Moreover, we write the result in terms of the 
physically relevant timescale in spin relaxation, r a , which 
for the case of interest here, u> <C T, is given by 



3m 



2tt 2 (1 - K) 



(44) 



In Ref. [3J] a single relaxation time approximation was 
made for the collision term, and this led to the response 
function (in the absence of mean-field effects) 



X 



relax 



N(0) 

1 — \LOT a 



In terms of t ct , the quantity £ is given by 



c 2 



(1 Ajx) 



(45) 



(46) 



and this ensures that the relaxation time approximation, 
Eq. (|4"5)) . reproduces the leading term \N(0)/(u)T a ) in the 
response function for lot„ 3> 1. This can be shown by ob- 
serving that, in this limit, the leading contribution to the 
deviation function is given by £ {0) (x) = l/[2cosh(x/2)]. 
When ujtq is large, the leading contributions to due to 



collisions can be calculated in perturbation theory, treat- 
ing l/(u>To) as small parameter, and one finds 



2ujtq 
A 



(xl+^)^( Xl ) 

dxidx-^dx^ 8 X 
_ OQ 4cosh(x2/2) cosh(x4/2) 



e (0) (x 3 ) 



2lutq 



(x 2 +7r 2 )(l-A ff )e (0) (^i), (47) 

where we have used the fact that l/[2cosh(x/2)] is the 
eigenfunction of the integral operator with eigenvalue 1, 
see Eq. (jAlOl) . From this, one finds that leading con- 
tribution to the response function is given by X a — 
\N(Q)/(uiTa) at high frequency. The analogous result for 
spin-independent distortions of the Fermi surface was de- 
rived in Ref. HI. 



IV. EXAMPLES 

One case that can be solved simply is that for X a = 0, 
and Eq. (fTS"]) then gives 



fa) 



1UT 



(x 2 + 7r 2 C 2 )cosh(a;/2) 



(48) 
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FIG. 1: (Color online) Left panels: Real and imaginary parts of the response function X„ in units of the density of states at 
the Fermi surface N(0), for \ a = 1/2,0 and —1/3, as a function of ujr a . The spin relaxation time r CT is given by Eq. (|44[) . 
For comparison, we show the response function based on the relaxation time approximation XJ elax . Right panels: Relative 
deviations of the real and imaginary parts of the relaxation time approximation and the exact result, (ReXJ elax /ReX CT ) — 1 and 
(ImXJ ela5 7lmX CT ) - 1, for the same values of \ a . 



Substituting this expression into Eq. 
response function [20 ] 



N(0) 1 



7T 2 C 



^'((1 + 0/2) 



301) leads to the 



(49) 



where tp'i 2 ) is the trigamma function [211 ]. In this limit, 
numerical results using Eq. (|49[) agree very well with 
those obtained from the general expression, Eq. (|4"3")l . 
thereby providing a check on the latter. 

In the middle- left panel of Fig.[TJ we show the real and 
imaginary parts of the response function X„ for \ a = 
as a function of ujt„. For comparison, results based on 
the relaxation time approximation, Eq. (|45[) . are also 
given. The physically relevant time scale that determines 
the transition between collisionless and hydrodynamic 
behavior is t„ which, according to Eq. (|4"4"|) , is shorter 
than To by a factor 3/(27r 2 ) « 0.15 for X a — 0. In the 
right panel of Fig. [TJ the relative differences between the 
relaxation time approximation and the exact result are 
shown more clearly, by plotting (ReX^ lax /ReX CT )-l and 
(ImA^ olax /hx\X a ) — 1. In the top and bottom panels of 
Fig. [TJ we present results for two other values of A CT , rel- 
evant for neutrino interactions and the spin response of 
neutron matter. These are X a = 1/2, which corresponds 
to a ratio of the quasiparticle relaxation time to the spin 
relaxation time, t/t ct = 3ro/(27r 2 r (T ) = 1/2, typical for 



low-momentum two-nucleon interactions (see Fig. 3 in 
Ref. 3) with very similar rates obtained in chiral effec- 
tive field theory [22[ , and A,j = —1/3, which corresponds 
to T/T a = 4/3, the value for the one-pion exchange ap- 
proximation to nuclear interactions Q ■ 

At low frequencies the relaxation time approximation 
always underestimates ImA CT . This may be understood 
in terms of the standard variational calculation of the 
relaxation time in the hydrodynamic limit. Generally, 
one can write the response function at low frequencies as 



X„ 



iV(0)(l 



hydro \ 



(50) 



where r^ ydro is a characteristic relaxation time in this 
limit. According to Eqs. (|4^|) and (|5TJ|) , the ratio of the 
spin relaxation times in the hydrodynamic and collision- 
less limits is given by 

T hydro _ 4(1 _ A ^ 



E 

r=0 



4r + 3 



(2r + 1) (2r + 2) [(2r + l)(2r + 2) - 2A fl 



(51) 



This result may also be obtained directly by generalizing 
the calculations of Refs. [H, Q to the problem of spin 
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FIG. 2: Ratio of spin relaxation times in the hydro dynamic 
and collisionless limits T% ydr ° /r a as a function of A CT . 

relaxation in a system with noncentral interactions. The 
quantity T a is the simplest variational estimate of the 
relaxation time and consequently always underestimates 
the relaxation time r^ ydro in this limit. In Fig. [21 we show 
the ratio r^ ydro /r cr as a function of \ a . For systems, in 
which collisions only occur between particles with the 
same spin projection and always change the sign of the 
spin projection, one has A CT = — 1 and the lack of spin 
conservation is maximal. In this case, T (T hydro /T (T « 1.17, 
while for A CT = the ratio is r CT hydro /T CT = tt 2 /9 ~ 1.10 
and for A CT -> 1 one finds r CT hydro /r CT -> 1. 

The relaxation time approximation overestimates 
lmX a at higher frequencies. That the differences of the 
imaginary parts cannot have the same sign for low and 
high frequencies is a necessary consequence of the fact 
that the response functions satisfy a Kramers-Kronig re- 
lation, and therefore 




= RcX a (0) - ReA^ olax (0) = . (52) 

In addition, at high frequencies KeX a is always greater 
than the relaxation time approximation. This is due to 
the fact that the leading term at high frequencies is pro- 
portional to the average of the square of the collision op- 
erator, while in the relaxation time approximation this is 
replaced by the square of the average. For X a = 0, this 
ratio can be calculated simply, and one finds 

Re ^ci ax (n<M* 2 +- 2 )sech 2 (./ 2 )y 5 

ReX a * /!° dx{x 2 + 7r 2 ) 2 sech 2 (a;/2) ~ 6 ' 

(53) 

in the limit wr CT 3> 1. 

In summary, our results show that the differences be- 
tween the real and imaginary parts of the spin response 
function calculated in the relaxation time approximation 
and the exact result are generally less than 10% for the 



conditions relevant to neutrino interactions in stellar col- 
lapse, supernovae and neutron stars. We do not expect 
this conclusion to change when mean-field effects are in- 
cluded, since for the densities of interest, the Landau 
parameter Go « 0.8 is positive and not large (see, for 
example, Ref. [23]). Consequently, the errors of using 
the relaxation time approximation are small compared 
to the theoretical uncertainties in determining the spin 
relaxation times in dense matter from nuclear forces. 



V. DISCUSSION 

We first consider applications of the spin response 
function to calculations of neutrino-pair bremsstrahlung 
and absorption and to neutrino scattering. As shown 
in Ref. [2J], if one neglects the momentum transfer to 
the nuclear medium and Pauli blocking effects for neutri- 
nos, the rates of neutrino processes when averaged over 
a Boltzmann distribution for the neutrinos are propor- 
tional to integrals of the form 

I»=] o duu 1 _ e _ u/T , (54) 

where the exponent ranges from v = for neutrino 
scattering to v = 7 for the rate at which neutrino- 
pair emission and absorption equilibrate the temperature 
of the nucleons with the temperature of the neutrinos 
and antineutrinos. For v = and provided the width 
(1 + Go)/t ct is small compared with T (which is a neces- 
sary condition for the picture of long-lived quasiparticle 
excitations to be applicable, see also Ref. Q), one can 
expand the denominator to first order in w/T and the 
rate for neutrino scattering is then proportional to 

r^^M {Q) (55) 

Jo u 

where the latter expression is a consequence of the 
Kramers-Kronig relation. Therefore, for this case the dif- 
ferences between the frequency dependences of xjj elax and 
Xa- have no effect. For larger values of v, higher frequen- 
cies are weighted more heavily than lower ones, and con- 
sequently the relaxation time approximation will overes- 
timate rates. However, the error can never be greater 
than the maximum deviation of Ivax™ from Imxa- 

There are a number of directions in which this work 
could be extended. One is to consider the case when 
u> > T, and another is to treat partially degenerate or 
nondegenerate nucleons. 

VI. CONCLUDING REMARKS 

In this paper, we have solved the quasiparticle trans- 
port equation for the spin response of a normal Fermi 
liquid for arbitrary values of u>T a . In a field-theoretical 
framework, this amounts to solving the Bethe-Salpeter 
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equation for a quasiparticle-quasihole pair excitation, 
taking into account self-energy and vertex corrections at 
the level of two-body collisions. One of the advantages 
of the transport equation approach is that conservation 
laws are properly taken into account. 

In particular, we have studied the relaxation of the 
isotropic part of the quasiparticle distribution function 
in momentum space. However, the basic formalism can 
be applied to relaxation of other spherical harmonic com- 
ponents, with appropriate changes in the definition of the 
quantity corresponding to \ a . These methods may also 
be applied to multicomponent systems, the main differ- 
ence being that the coefficients a r then depend on the 
species, and one must solve a set of linear equations. The 
techniques therefore represent a generalization of those 
used to calculate transport coefficients of multicompo- 
nent Fermi liquids [25j ]. 

The differences between the exact solution and the re- 
laxation time approximation are small compared with the 
uncertainty of other input to calculations of neutrino pro- 
cesses in dense matter. However, the effects could be 
of importance in related situations, such as the absorp- 
tion of sound in liquid 3 He for conditions between the 
collisionless and hydrodynamic regimes. Another poten- 
tial application is to calculate the frequency-dependent 
conductivity of a charged Fermi liquid and the magne- 
toresistance of metals, problems attacked in Ref. p} by 
variational techniques. 
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APPENDIX A: MATHEMATICAL DETAILS 

To make the paper reasonably self-contained, we list 
a number of useful mathematical relations in this ap- 
pendix. In the derivation of the quasiparticle transport 
equation, Eq. (fT8|) . we encountered the integrals 

K (x) - [°° dx2dxi S ( x + X2 ~ x ^ (All 
~~ 7-oc 4cosh(iE2/2)cosh(iE4/2) ' ( ' 



and 



K 3 (x) = 



dxidxzdx^ S(x + X2 — x 3 — X4) 
cosh(x2/2) cosh(x3/2) cosh(a;4/2) 



These may be evaluated by representing the delta func- 
tion as 

dt - t „ 



S(x) 



2tt 



(A3) 



and using the methods described in Appendix A of Chap- 
ter 1 of Ref. UJl . The results are 



K 2 (x) = 



2sinh(:r/2) ' 



and 



K 3 (x) = 



2cosh(cc/2) 



(A4) 



(A5) 



By multiplying Eq. (|3"4"|) by <Ev<(fc), subtracting the 
same expression with r and r' interchanged, and then 
integrating over k, one sees that the $ r (/c) are orthogo- 
nal in the sense that 



$ r (AQ gy(fc) _ U r 

Qif\> Off 1 ■ 

cosh (irk) 27r 



(A6) 



where U r is given explicitly in Eq. (|41[) . For the eigen- 
functions in terms of x, this condition is equivalent to 



dx {x 2 + 7r 2 C 2 ) ® r (x) $ r '(x) = 8 rr < 



U r A r 



(A7) 



(A2) 



In the hydrodynamic limit (£ — > 1), the eigenfunctions 
Q r (x) have a simple form, since 

r°° dk ~ 
®r{x)= $ r (fc)e ifc *, 

J -oo 27T 

/OO T 1 

sech(^fe) CU 2 (tanh(Trfc)) e ikx , (A8) 

where we have used Eq. (|37| for $ r (fc). The integral can 
be evaluated by contour integration, closing the contour 
in the upper half plane for x > and in the lower half 
plane for x < 0. The integrand has poles of order 2r+l at 
k = i(^+l/2) with integer v, and therefore the integral is 
proportional to the sum of the residues of the integrand 
at these points, which is a geometrical progression with 
the ratio of subsequent terms equal to — e - '^'. The final 
result is $ r (x) ~ Q r (x)/[2 cosh(x/2)], where Q r (x) is 
the residue of e lkx C 2 i / 2 (coth(7rfc))/ sinh(7rfc) at k = 0, a 
polynomial of order 2r. The polynomials Q r {x) satisfy 
the orthogonality relation 

x 2 + 7T 2 

dx 7i Qr(x) Q r i(x) = , for r^r'. 

cosh 2 (.T/2) V ; w 

(A9) 

(A10) 
(All) 

(A12) 



The first few eigenfunctions are given by 
Qo(x) - 1 , 



Qi(x) 
Q 2 (x) 



1-^). and 



1 - 



14 x 2 
~5 V* 



15 V 4 
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For other values of C the eigenfunctions do not have such a simple form. 
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